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I. ABSTRACT. 

Numerical results based on an extended BPM 
algorithm indicate that, in Marcatili's lossless ta- 
pers and bends, through- flowing waves are dras- 
tically different from standing waves. The source 
of this surprising behavior is inherent in Maxwell's 
equations. Indeed, if the magnetic field is correctly 
derived from the electric one, and the Poynting 
vector is calculated, then the analytical results are 
reconciled with the numerical ones. Similar con- 
siderations are shown to apply to Gaussian beams 
in free space. 



II. INTRODUCTION. 

In 1985, Marcatili jlj] infringed a historical 
taboo, by showing that lossless tapers and bends 
in dielectric waveguides can be conceived and de- 
signed, at least on paper. The key feature shared 
by all the infinity of structures which obey Mar- 
catili's recipe, is the fact that the phase fronts of 
the guided modes which propagate in them, are 
closed surfaces. As well known, phase fronts which 
extend to infinity in one direction orthogonal to 
that of propagation do entail radiation loss, but 
closed fronts can avoid this problem. However, 
shortly after the first recipe M , it was pointed out 
H that that recipe could generate some inconsis- 
tencies. In fact, a traveling wave with a closed 
phase front is either exploding from a point (or a 
line, or a surface), or collapsing into such a set. 
In a lossless medium where there are no sources, 
this is untenable. On the other hand, it was also 
pointed out in Q that a standing wave with closed 
constant-amplitude surfaces is physically meaning- 
ful. Therefore, propagation of a through-flowing 
wave through any of Marcatili's lossless tapers or 
bends has to be described in this way: the incom- 
ing wave must be decomposed as the sum of two 
standing waves, of opposite parity with respect to a 
suitable symmetry surface. The output wave was 
then to be found as the sum of the values taken 
by the two standing waves at the other end of the 



device. Another point raised in ^ was that very 
similar remarks apply to Gaussian beams in free 
space. 

Later on, the literature showed that interest in 
this problem was not so high, for a long time. Re- 
cently, though, we observed several symptoms of a 
renewed interest in low-loss and lossless |].[7j 
tapers or bends. This induced us to try to go be- 
yond the results of and to clarify further the 
difference between through-flowing and standing 
waves in Marcatili's tapers. 

The new results reported in this paper can be 
summarized as follows. In Section III, we show 
that the numerical analysis (based on an extended 
BPM algorithm) of Marcatili's tapers reconfirms 
that indeed through-flowing waves are drastically 
different from standing ones. The latter ones 
match very well the analytical predictions of the 
original recipe , but through-flowing waves have 
open wave fronts, which do not entail any phys- 
ical paradox. In Section IV, we provide an ana- 
lytical discussion of why, in contrast to what oc- 
curs with plane waves in a homogeneous medium 
and with guided modes in longitudinally uniform 
waveguides, through-flowing waves are so differ- 
ent from standing ones. We show that this is a 
rather straightforward consequence of Maxwell's 
equations. From this we will draw the conclu- 
sion that a through-flowing wave propagating in 
one of Marcatili's tapers is never strictly lossless. 
Nonetheless, our numerical results reconfirm that 
the recipes given in [Q do yield extremely low ra- 
diation losses. 

Finally, we address briefly the case of Gaussian 
beams in free space, and explain why they be- 
have essentially in the same way as the devices 
we discussed above. In fact, Maxwell's equations 
show that in general the phase fronts of the mag- 
netic field in a Gaussian beam are not the same 
as the phase fronts of the electric field. Therefore, 
the Poynting vector is not trivially proportional to 
the square of the electric field. Consequently, a 
through-flowing beam, resulting from two super- 
imposed standing waves of opposite parities, can 
be surprisingly different from the parent waves. 



III. NUMERICAL RESULTS. 

The geometry of Marcatili's tapers can eventu- 
ally be very complicated (e.g., see ||). For our 
tests, however, we chose a simple shape, to avoid 
the danger that geometrical features could hide the 
basic physics we were trying to clarify. The re- 
sults reported here refer to a single-mode taper 
whose graded-index core region is delimited by the 
two branches of a hyperbola (labeled A and A' in 
Figs. 1 and 2), and has a mirror symmetry with 
respect to its waist. This is a "superlinear" taper, 
according to the terminology of |jlj, with an index 
distribution (see again B) 
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where rj and $ are the elliptical coordinates, in 
the plane of Figs. 1 and 2. Fig. 1 refers to a stand- 
ing wave of even symmetry with respect to the 
waist plane, Fig. 2 to a standing wave of odd sym- 
metry. The closed lines are constant-amplitude 
plots. We see that they are essentially elliptical, 
so they agree very well with the predictions of M. 




FIG. 1. Constant-amplitude plot of a standing wave 
of even symmetry (with respect to the waist plane) in 
a superlinear Marcatili's taper. 

As mentioned briefly in the Introduction, these 
results were generated using an extended BPM, 
which deserves a short description. In fact, it is 
well known that the standard BPM codes are suit- 
able to track only traveling waves, as they neglect 
backward waves. Our code (using a Pade's op- 
erator of order (5,5)) also generates a traveling 
wave, and its direction of propagation is inverted 
whenever the wave reaches one of the taper ends. 
In order to generate single-mode standing waves, 
each reflection should take place on a surface whose 
shape matches exactly that of the wave front. This 



is very difficult to implement numerically, but the 
problem can be circumvented, letting each reflec- 
tion take place on a phase-conjugation flat mirror. 
Our code calculates then, at each point in the ta- 
per, the sum of the forward and backward fields, 
and stops when the difference between two itera- 
tions is below a given threshold. 




FIG. 2. Constant-amplitude plot of a standing wave 
of odd symmetry (with respect to the waist plane) in 
a superlinear Marcatili's taper. 

Figs. 3 and 4 refer to a through-flowing wave. 
The almost horizontal dark lines in Fig. 3 are its 
phase fronts. They are drastically different from 
those predicted by the analytical theory in |5J], 
which are exemplified in the same figure as a set 
of confocal ellipses. Note that the through-flowing 
wave has been studied numerically in two ways. 
One was simply to launch a suitable transverse 
field distribution, and track it down the taper with 
a standard BPM code. The other one was to cal- 
culate the linear combination (with coefficients 1 
and j) of the even and odd standing waves shown 
in Figs. 1 and 2. The results obtained in these two 
ways were indistinguishable one from the other. 
This proves that indeed through-flowing waves are 
drastically different from standing ones. In par- 
ticular, as we said in the Introduction, they keep 
clear from any paradox connected with energy con- 
servation. 

Fig. 4 shows a field amplitude contour plot for 
the same through-flowing wave as in Fig. 3. It 
indicates that propagation through the taper is in- 
deed almost adiabatic. Therefore, as anticipated 
in the Introduction, insertion losses of Marcatili's 
tapers are very low (at least as long as the length 
to width ratio is not too small), although they are 
not strictly zero. As a typical example, for a total 
taper length of 2.5/im, a waist width of 0.55/im and 
an initial-final width of 1.65/im, BPM calculations 
yield that the lost power fraction is 1.4 x 10~ 4 . A 
typical plot of power vs. distance along a taper 



with these features is shown in Fig. 5. 
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FIG. 3. Phase fronts of a through-flowing wave in 
the same superlinear taper as in Figs. 1 and 2. 

IV. THEORETICAL DISCUSSION. 

For the sake of clarity, let us restrict ourselves to 
the case of two-dimensional tapers, like those of the 
previous section, where the geometry and the index 
distribution are independent of the z coordinate, 
orthogonal to the plane of the figures. However, 
our conclusions will apply to 3-D structure also. 

The index distributions found in the corner- 
stone paper jl] are such that the TE modes (elec- 
tric field parallel to z) satisfy rigorously a wave 
equation which can be solved by separation of vari- 
ables. Obviously, the same equation is satisfied rig- 
orously by the transverse component of the mag- 
netic field. 

However, in general two solutions of these two 
wave equations which are identical, except for a 
proportionality constant, do not satisfy Maxwell's 
equations in that structure. This is very easy to 
show, for example, for the case which was called 
"linear taper" in m, namely, a wedged-shape re- 
gion with a suitable index distribution, where a 
guided mode propagates in the radial direction. 
The claim [0 that the dependence of E z on the 
radial coordinate is expressed by a Hankel func- 
tion of imaginary order iv, related to other fea- 
tures of the taper, is perfectly legitimate. How- 
ever, one cannot extrapolate from it that the same 
is true for the magnetic field. In fact, calculat- 
ing the curl of the electric field we find that the 
azimuthal component of the magnetic field is pro- 
portional to the first derivative of the Hankel func- 
tion, which is never proportional to the function 
itself. The same is true for the Mathieu function 
of the fourth kind, which satisfy the wave equation 
in the coordinate system which fits the superlinear 
taper of the previous Section. This entails a dras- 



tic difference with plane waves, and with guided 
modes in uniform waveguides, where the deriva- 
tive of the exponential function that describes the 
propagation of the electric field is proportional to 
the function itself. In the cases at hand, the con- 
cept of wave impedance becomes ill-grounded. In 
fact, the electric field and the transverse magnetic 
field have identical dependencies on the transverse 
coordinate, so that their ratio is constant over each 
wavefront, but they are different functions of the 
longitudinal coordinate, as if the 'wave impedance' 
were not constant at all along the wave path. This 
indicates why it is very risky, in the case at hand, to 
make general claims on the Poynting vector start- 
ing from the spatial distribution of only the electric 
field. To strengthen our point, let us prove explic- 
itly that it is not self-consistent to claim that a 
purely traveling TE wave, whose radial dependence 
is expressed by a Hankel function of imaginary or- 
der, Hi Vl can propagate along a linear taper. As 
we just said, for such a wave E z is proportional 
to H iv , is proportional to H' iv , so the radial 
component of the Poynting vector is proportional 
to Hi„(H' iv )* . In a purely traveling wave there is 
no reactive power in the direction of propagation. 
Combining with what we just said, it is easy to see 
that this would imply \H^\ 2 — constant along the 
(radial) direction of propagation, a requirement 
that is not satisfied by Hankel functions. (Note, 
once more, that it is satisfied by exponential func- 
tions). Therefore, any wave along a linear taper 
whose radial dependence is expressed as a Hankel 
function must be at least a partially standing wave. 
A through-flowing wave, if it exists, must behave 
in a different way. 
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FIG. 4. Field-amplitude contour plot, in the same 
superlinear taper as in the previous figures. 

Finally, let us address briefly the case of Gaus- 
sian beams in free space. It was pointed out in 
g that they behave essentially in the same way as 
the devices we discussed above. There is still some- 
thing to add to the discussion of M . Assume that 



the electric field of an electromagnetic wave has the 
classical features of a TEM 00 Gaussian beam (see, 
e.g., ||). Then, Maxwell's equations show that the 
phase fronts of the magnetic field are not the same 
as those of the electric field, neither on the waist 
plane nor far from it. Hence, the Poynting vector is 
not trivially proportional to the square of the elec- 
tric field. This entails the presence of a reactive 
power (never accounted for in the classical class- 
room explanations of Gaussian beams), and an ac- 
tive power flow which is not always along the lines 
orthogonal to the electric field phase fronts. Once 
again, a through-flowing beam, resulting from two 
superimposed standing waves of opposite parities, 
is different from the parent waves, and the differ- 
ence is maximum on the symmetry plane, i.e. at 
the beam waist. Due to time and space limits, the 
details of this discussion must be left out of this 
presentation, and will be published elsewhere j[o| . 




FIG. 5. Power vs. distance, in a superlinear taper 
of the shape shown in the previous figures, whose pa- 
rameters are specified in the text. 

V. CONCLUSION. 

We tried to shed new light on an old problem, 
namely, whether the idea of a guided mode trav- 
eling without any loss through a dielectric taper 
can be sustained without running into any physical 
paradox. Our numerical results, obtained with an 
extended BPM technique, have fully reconfirmed 
what was stated in Q: in Marcatili's tapers, stand- 
ing waves have the basic features outlined in [jjj, 
but through-flowing waves do not. This prevents 
them from running into a paradox, but on the 
other hand entails some loss, although very small 
indeed. Next, we have provided an explanation 
for the unexpected and puzzling result, a drastic 
difference between standing and through-flowing 
waves in the same structures. The source of these 
"surprise" is within Maxwell's equations. 

It was pointed out in [Q that some of the prob- 
lems discussed here with reference to Marcatili's 



tapers apply to Gaussian beams in free space as 
well. 
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